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O ■ ABSTRACT 

oo ■ 
O 

The COMPLETE Survey of Star-Forming regions offers an unusually comprehensive 
)«j ■ and diverse set of measurements of the distribution and temperature of dust and gas in 

^ ■ molecular clouds, and in this paper we use those data to find new calibrations of the 

^ ■ "X-factor" and the ^^CO abundance within Perseus. To carry out our analysis, we: 1) 

apply the NICER (Near-Infrared Color Excess Method Revisited) algorithm to 2MASS 
data to measure dust extinction; 2) use dust temperatures derived from re-processed 
IRAS data; and 3) make use of the ^^CO and ^^CO (1-0) transition maps gathered by 
COMPLETE to measure gas distribution and temperature. Here, we divide Perseus 
into six sub-regions, using groupings in a plot of dust temperature as a function of 
LSR velocity. The standard X factor, X = N{Yl2) /W{^'^CO), is derived both for the 
whole Perseus Complex and for each of the six sub-regions with values consistent with 
previous estimates. However, the X factor is heavily affected by the saturation of the 
emission above Ay^4: mag, and variations are found between regions. We derive linear 
fits to relate W{^'^CO) and Ay using only points below 4 mag of extinction. This 
linear fit yields a better estimation of the Ay than the X factor. We derive linear 
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relations of W(^^CO), N(^^CO) and W{C^^O) with Ay- In general, the extinction 
threshold above which ^'^CO (1-0) and C^^O (1-0) are detected is about 1 mag larger 
than previous estimates, so that a more efficient shielding is needed for the formation 
of CO than previously thought. The fractional abundances (w.r.t. H2 molecules) are 
in agreement with previous works. The (1-0) lines of ^^CO and ^^CO saturate above 
4 and 5 mag, respectively, whereas C^^^O (1-0) never saturates in the whole Ay range 
probed by our study (up to 10 mag). Approximately 60% of the positions with ^^CO (1- 
0) emission have sub-thermally excited lines, and almost all positions have ^'^CO (1-0) 
excitation temperatures below the dust temperature. We compare our data with PDR 
models using the Meudon code, finding that ^^CO (1-0) and ^^CO (1-0) emission can 
be explained by these uniform slab models with densities ranging between about 10^ 
and 10^ cm~'^. In general, local variations in the volume density and non-thermal 
motions (linked to different star formation activity) can explain the observations. Higher 
densities are needed to reproduce CO data toward active star forming sites, such as 
NGC 1333, where the larger internal motions driven by the young protostars allow more 
photons from the embedded high density cores to escape the cloud. In the most quiescent 
region, B5, the ^^CO and ^^CO emission appears to arise from an almost uniform thin 
layer of molecular material at densities around 10^ cm~^, and the integrated intensities 
of the two CO isotopologues are the lowest in the whole complex. 

Subject headings: dust, extinction — ISM:abundances — ISM:molecules — ISM:individual 
(Perseus molecular complex) 



Introduction 



Although H2 is the most abundant molecule in the interstellar medium (by about four orders 
of magnitude), it cannot be used as a tracer of the physical conditions in a molecular cloud. 
In fact, being a homonuclear species, II2 does not have an electric dipole moment and even the 
lowest (electric quadrupole) rotational transitions require temperatures and densities well above 
those found in typical mol e cular clouds. The next most abundant molecule is ^^CO, and since its 



discovery by lWilson et a. 



(I1970I) it ha s been considered the best tracer of H2 and the total mass 



of molecular clouds (e.g. ICombeslll99ll ). Because of its relatively large abundance and excitation 
properties, the ^^CO (1-0) line is typically optically thick in molecular clouds, so that rarer CO 
isotopomers (in particular ^^CO) have also been used to trace the cloud mass in the most opaque 
regions. However, previous attempts to derive cloud masses from the thinner isotopomers present 
a large scatter (factor of 10), indicating that ^^CO line strength is not an entirely straightforward 
indicator of gas column. In particular, the following points should be taken into account to derive 
less uncertain conversions: 1) preferential photodissociation of ^^CO at low optical extinctions. Ay] 
2) active chemical fr actioiiation in the presence of ^'^C^ ions, enhancing the relative abundance of 
^^CO in cold gas ( WatsonI 119771 ): and 3) variations in the ^^CO/^^CO abundance ratio across the 
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Galaxy. 

The conversion between ^^CO (1-0) integrated intensity {W(^'^CO)) and H2 column density 
(A^(H2)) is usually made using the so-called "X-factor," defined as 



X = 



iV(H2 



vF(i2co) ■ 



(1) 



In order to calibrate this ratio, one needs to measure the column of H2 which is difficult to do 
reliably. 

One method to derive this ratio use s only the ^•^CQ e missi on and the assumption that molecular 
clouds are close to virial equilibrium. I Solomon et al.l (|l987l ) found a tight relation between the 
molecular cloud virial mass and the ^^CO luminosity (Myx = 39(Lco)'^''^^)- This relation enabled 
them to derive a ratio of 3.0 x lO^^^ cm^^ K^^ km^^ s for the median mass of the sample (10^ Mq). 
In this case, the main source of uncertainty is the assumption of virial equilibrium for the molecular 
clouds. 

Another method uses the generation of gamma rays by the collision of cosmic rays with H 
and H2. The H2 abundance is calculated us ing (assumed opticall y thin) 21-cm observations of 
H I and gamma ray observations in concert. iBloemen et al.l (119861) using COS-B data d erived a 
value of (2.8 ±0.4) x lO^o cm-^ K'^ km"^ s for the X factor, while Istrong Matto^j (|l996l ) derived 



(1.9 lb 0.2) X 10^'^ cm^^K^^km^^s using the EGRET data. Uncertainties in gamma-ray-derived 
estimates of the X-factor stem from coarse resolution (0.5° and 0.2° for COS-B and EGRET, 
respectively) and from the assumption that there are no point sources of gamma rays. 

An alternative method to derive the column density of H2 makes use of the IRAS 100 fim data 
to estimate the total column density of dust, which can then be used, assuming a constant dust- 
to-gas ratio (=0.01, value adopted also here), to estimate the total column density of Hydrogen 
(N£H) = iV(H I) -|- 2A^(Il2)) and derive the conversion factor X. Using this method iDame et al. 
(l200lh derived a valu e of (1.8 ± 0.3) x 10^° cm^^ j^-i g ^j. ^j^g jjg]^ ^.^e Milky Way, while 



de Vries et al. (jl987l ) derived ( 0.5 ifc 0.3) x 10^ ° cm~2 km~^ s for the high-latitude far- infrared 
"cirrus" clouds in Ursa Major. iFrerking et al.l ( 19821 ) find no correlation at all between W{^'^CO) 
and A'^(H2) in Taurus (where the ^^CO (1-0) integrated intensity is roughly constant above 2 mag 
of visual extinction) while in p-Ophiuchus they derive (1.8 it 0.1) x lO^'^ cm"^ km" 



s. 



More recently, iLombardi et al.l (|2006l ) studied the Pipe Cloud using an extinction map derived 
using the Near Infrared Color Excess Revisited (NICER) technique on 2MASS and ^^CO (1-0) 
data. They derived an X factor of (2.91 ± 0.05) x 10 



2° cm 2 



K ^ km ^ s, similar to what is found 



in cloud core reg ions and dark nebulae (see compilation by lYoung &: Scovillelll982l ). However, the 
factor derived by lLombardi et al.l (j2006l ) includes a correction for Helium and uses the non-standard 
expression A^(H) = iV(H l)+N{R2), while using the standard definition of iV(H) = A^(H l)+2N(H2) 
they would derive (1.06 ± 0.02) x lO^o cm"^ j^-i j^j^- - 



s. 



Frerking et al.l ( 19821 ) studied the relation between visual extinction and ^^CO data, and found 
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that the number of H2 molecules per i^CO molecule, [Ha/^^CO], is 3.7 x 10^ in Taurus and 3.5 x 10^ 
in Ophiuchus, with extinction thresh olds (below wh i ch no ^^CO is detected) of 1.0 and 1.6 mag, 
respectively. To make these estimates, iFrerking et al.l (|l982l ) derived pencil-beam extinctions along 
the line of sight towards a handful o f positions with background stars, by using near- infrared 
(NIR) spectroscopy from lEliad (|l978l ) to estimate the spectral type and derive the extinction. 
However, pencil beam observations do not trace the same material pro bed by the molecular line 
emission observations, and this can introduce large uncertainties. In fact. lArce &: GoodmanI ()1999l ) 
compared spectroscopically-determined extinction and IRAS-derived extinction in a stripe through 
Taurus, finding a 1-a dispersion of 15% when the extinctions are compared between 0.9 < Ay < 
3. mag. This unavoida ble dispersion is likely to affect previous column density derivations, such 



as 



Frerking et al.l (|l982l ). as well. 



Previous studies of the ^"^CO abundance have also been carr ied out specifically in the por - 
tions of the Perseus Molecular Cloud Complex we study here. iBachiller Cernichard (119861 ) 
used an extinction map deri ved from star counting on the Palomar Observatory Sky Survey 



(jCernicharo Sz Bachilled Il984l ). The angular resolution of their extinction map is 2.5', which is 
smoothed to 5' resolution for comparison with the 4.4' resolution of their ^^CO data, from which 
they d erive an [Ha/^ '^ COI abundance and threshold extinction of 3.8 x 10^ and 0.8 mag, respec- 
tively^ Laji2ei_et_jl. J ^89) compare ^^CO data of 2' resolution with the extinction map derived 
by ICernicharo k BachilleJ |l984) , obtaining [H2/^^C0]= 3.6 x 10^ and an extinction threshold of 
0.5 mag. However, this study was done only in the B5 region. A summary of earlier measurements 
of the X-factor and ^^CO abundances can be found in 16.21 and! 



The main disadvantage of previous calibrations done on Perseus is that t hey use optically-hased 
star counting to create ex tinction maps of regions with h i gh vis ual extinction (jBachiller &: Cernicharo 
1986l : lLanger et al.lll989l'). For Qphiuch us JSchnee et al.l (j2005l ) find that the extinction derived using 
optical star counting ( Cambresvlll999l ) is systematically underestimated by ~ 0.8 mag in compari- 



son with near-IR-based extinction mapping. The offset (and resulting inaccuracy) is caused by the 
difficulty in fixing the zero-point extinction level - in such a high extinction region - in the optical 
star counting extinction map. In addition, the derived Ay from optical star counting in Perseus 
does not have a large dynamic range {Ay < 5 mag), while the NICER extinction niap use d in our 
work is accurate up to 10 mag with small errors (cxyiv < 0-35 mag; see iRidge et al.ll2006al ). 



The B5 region in Perseus is special, in that it has been the subject of an unusually high number 
of studies seeking to understand its basic physical properties. As shown in Figure [3l B5 is somewhat 
isolated from the rest of Perseus, and it is o nly forming a very small number of stars, making it 
an obviously good choice for detailed study. lYoung et al.l (119821 ) used a LVG model to derive an 
average de nsity of 1.7 x 1 0^ cm~^ and kinetic temperatures ~ 10 — 15 K, but only for some stripes in 
the cloud. iBenschI ((20061) used ^^CO and ^^CO maps with C I pointing observations of 12 positions 
in a North-South stripe from the central B5 region to model the emission with a PDR code. From 
this analysis he derives average densities ~ 3 X 10^ - 3 X 10"^ cm . 
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The COMPLETE dataset offers the unique opportunity to study the emission of various CO 
isotopologues across the whole Perseus complex with unprecedented sensitivity and spatial resolu- 
tion. In the present paper, COMPLETE data are analysed in detail to measure the ^^CO excitation, 
the ^^CO abundance and the X factor across the complex and study their variations. The observed 
changes in the measured quantities are then related to local properties of the gas and dust. Using 
PDR codes, we find that local variations in the volume density and non-thermal motions (linked 
to different star formation activity) can explain the observations. 

The extinction map, molecular and IRAS data used in this paper are presented in Sect. [2j 
The data selection is discussed in Sect. [3l The six regions in which the Perseus Molecular Cloud 
Complex is divided are identified in Sect.Sl Section [5] contains the analysis of the data, including 
the ^^CO column density determination and the curve of growth. Results can be found in Sect. [HI 
The comparison between observations and PDR models is in Sec. [7| and conclusions are listed in 
Sect. El 



Data 



2.1. Extinction Map 



We use the Near-Infrared Color Excess Method Revisited (NICER) technique (jLombardi &: Alves 



200 ll ) on the Two Micron All Sky Survey (2MASS) point so urce catalog to calculate K-band extinc- 
tions. To derive Ay-, we use the relation Ak = 0.112 Ay (jRieke &: Lebofskvlll985l ). The resulting 
extinction map has a fixed r esolution of 5 ^ and the p ixel scale is 2.5^ The total size of the map is 
9° X 12° and is presented in lRidge et al.l (j2006al ) (see lAlves et al.l 120071 . for more details). 



Figure [U shows the estimated error (uncertainty) in the derived Ay for all the points in the 
extinction map used in this study. The uncertainties are correlated with the extinction values 
because more stars per pixel allow for a more accurate measure of extinction. Nevertheless, the 
linear slope in Figure [1] is smaller than 0.01, so that the fractional uncertainty in any pixel's 
extinction measure remains quite small compared to its value. The median of the error for all of 
Pe rseus using NICER on 2MASS data is 0.2 r nag, while in the case of the e x tinction map derived 
by ICernicharo &: Bachilleii (jl984l ) and used in iBachiller &: Cernichard ()1986l ) ; iLanger et al.l ()l989l ) 
the typical error is ^ 0.5 mag. 

While the improved accuracy offered by our new maps is significant, the increased dynamic 
range is even more critical to our analysis. The NICER extinction map o f Perseus is accurate up 
to Ay = 10 mag, while the extinction map derived using star counting by ICernicharo &: Bachiller 
(jl984l ) dies out above ~ 4 — 5 mag of visual extinction. 
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2.2. Molecular Data 



We use the line maps of the COMPLETE Survey (jRidee et alJl2006al l to estimate ^^CO (1-0) 
and ^^CO (1-0) column densities. Both lines were observed simultaneously using the FCRAO tele- 
scope. The line maps cover an area of ~ 6.25° x 3° with a 46" beam in a 23" grid, and the positions 
included in our analysis are shown in FigureO We correct the ^^CO (1-0) and ^^CO (1-0) maps for a 
main-beam efficiency assumed to be 0.45 and 0.49, respectively (http : //www, astro .umass . edu/$\sim$f cr ao/[ ). 
The flux calibration uncertainty is assumed to be 15% (Mark Heyer, private communication). 



In comparing the ^^CO integrated intensity presented in lDame et al.l ( 200ll ) with the integrated 
intensity of our data smoothed to match the 1/8° beam resolution, we find the COMPLETE 
data and Dame's measurements to be well fitted by a linear relation of slope 0.9 it 0.1 and off- 
set —2 lb 1 Kkr ns"^. We also com pare the COMPLETE ^^CO integrated intensity with data 
from Bell Labs (jPadoan et al.lll999l ). after smoothing our data to the 100" resolution and 1' grid 
of Bell Labs data. These data are fitted by a linear relation of slope 0.986 it 0.003 and off-set 
—0.53 lb 0.01 Kkms~^, where the small deviation from the 1 : 1 relation is most probably due to a 
small misalignment found between the images. 



In addition, we use the C^O (1-0) data-cube presented by Hatchell fc van der TakI (|2003l ) 



(and converted into FITS using CLASS90; iHilv-Blant et al. 
a smaller coverage and lower signal to noise. 



2005-1 



, taken with FCRAO but with 



We convolve all data-cubes with a Gaussian beam to obtain the same 5' resolution as the 
NICER map, and then re- grid them to the extinction map grid of 2.5'. 



2.3. Column Density and Dust Temperature from IRAS 



To estimate column density accurately from far-infrared (thermal) flux, one needs to also 
measure or calculate a temperature. Normally, this is accomplished by making measurements at 
two separated far-IR wavelengths, making assumptions about dust emissivity, and then calculating 
two "unknowns" (column density and temperature) from the two measurements of flux. In the 
happy case where extinction mapping offers an independent measure of column density over a wide 
region, the conversion of FIR flux ratios to column density and temperature can be optimized so as 
to minimize point-to-point differences in comparisons of extinction- and emission-derived column 
density. The column densities and temperatures derived f rom dust emission that we use in this 
paper and in lGoodman et al.l (|2007l ) come from the work of lSchnee et al.l (|2005l ). who re-calibrated 
IRAS-based maps by using the 2MASS/NICER extinction maps discussed above to constrain the 
column density conversions. 

hood) ssing 



The "IRAS" data used in 



Schnee et al 



of the IRAS Survey; iMiville-Deschenes &: Lagachd |2005| ) flux maps at 60 and 100 /im. The IRIS 
data have better zodiacal light subtraction, calibration and zero level determinations, and destriping 
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than the earlier ISSA IRAS survey release. 



There are two important caveats to apply to the FIR-based column density and temperature 
measurements we use here (applicable to previous work as well). First, the column densities are only 
optimized to reduce scatter in a global (Perseus-wide) comparison of dust extinction and emission 
measures of column density - they are still calculated based on the measured FIR fluxes at each 
point, and are thus not constrained to be identical to the 2MASS/NICER values. Second, the dust 
temperatures derived by this method are uncertain due to: 1) unavoidable line-of-sight temperature 
variations, which cause increased scatter in comparison s with extinction-ba sed measures and also 
cause a bias toward slightly higher temperatures (see Schnee et al.l 2006 ): and 2) the effect of 
emission from transiently heated Very Smal l Grains (seelSchnee et al.ll2008l ') . That said, the column 



density and temperature estimates based on lSchnee et al. 
a dramatic improvement. 



(|200d ) 



and used in this paper do represent 



3. Data Editing for Analysis 

The total number of pixels with data in all maps is 3765. But in order to have high-quality 
data in every pixel of every map used in our analysis, we trim our maps to exclude particular 
positions where any data are not reliable. The procedure used in data editing is described in this 
Section. 



3.1. Extinction Map 

Regions with both high stellar density and high extinction are typically associated with embed- 
ded populations of young stellar objects (YSOs). Thus, in creating extinction maps, including these 
regions, it would be foolish to assume that stars are background to the cloud and have "typical" 
near-IR colors. Among YSOs, the more evolved objects observable in J, H, and K bands could in 
principle still be used to measure the extinction, but the fact that they are not background objects 
still produces an underestimation of the extinction. This affects the precision of the method much 
more than the non-stellar IR colors (infrared excess) of YSOs. Therefore, we exclude all the regions 
with a stellar density larger than 10 stars per pixel from our analysis. This criterion removes 34 
pixels from our maps. In addition to this editing, we exclude a box around each of the two main 
clusters in Perseus: IC348 and NGC1333, removing another 158 pixels from our data, as is evident 
in Figure El (see Tabled]). 
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3.2. Molecular Transitions 

To remain in tlie analysis, lines must have positive integrated intensities in both ^^CO and 
^^CO and peak brightness temperatures of at least 10 and 5 times the RMS noise for ^^CO and 
^^CO respectively. 

Since ^■^CO is more abundant than ^'^CO, and self-shielding is more effective for the ^■^CO (1-0) 
transition than for the ^'^CO (1-0) transition, ^■^CO emission is always more extended than ^^CO, 
both spatially and kinematically. In addition, ^^CO lines are more affected (broadened) by outflows 
and are optically-thicker than ^^CO. Thus, line-widths for ^^CO, cr(^^CO), should always be larger 
than those of ^^CO, o"(^^CO). As a result, we keep only positions with 

a(i2co) > O.SaC^CO) , 

where the 0.8 factor has been chosen to take into account the uncertainties in the line- width 
determination. The line-widths and central velocities of the spectra are obtained through Gaussian 
fits. This filtering accounts for only 72 out of the 324 pixels edited out in this study. 

3.3. Final Data Set 

The result of our editing leaves us with 3400 pixels (shown in Figure [3l see § H] for details on 
color coding), out of an original 3724 pixels with ^^CO detections. Note that our pixels are 2x 
oversampled (see § 12.21) . so that 3400 pixels amounts to 3400/4 independent measures. 

4. Region Identification 

In trying to study Perseus as one object, it became apparent that much of the scatter in 
both the X-factor and in ^^CO abundance is caused by region-to-region variations (see §§ 16.21 and 
16. 3p . So, we have divided the Perseus Molecular Cloud Complex into six regions with the help 
of several plots comparing different parameters (e.g. ^^CO line-width, ^^CO LSR velocity, dust 
and excitation temperature). Figure [2] is an example of such plots, showing the ^^CO velocity, 
Vls'h(^^CO) (the Vlsb. of ^^CO and ^^CO are very similar) as a function of dust temperature, 
Td, for all the Perseus data. In this way, we avoid more arbitrary choices and minimize the data 
overlap among the different regions. As can be seen in Figure [21 the six regions cluster around 
characteristic values of Vlsr and T^, so they are easily identified. Minor further refinements on 
the region definitions is done to keep the regions physically connected, as shown in Figure [31 This 
step is needed primarily because there are regions with two velocity components along the line of 
sight and the Gaussian fitting can spontaneously switch from one component to the other. The 
effect of the two components is clearly seen in the points below 2.5 kms~^ in Figure[2l where points 
from three different geographical regions are merged into a single region of Vlsr ~ space. In 



FigureOwe show the final defined regions: B5, IC348, Sheh, Bl, N GC1333 and Westend . The Shell 
region is essentially the same as the shell-like feature discussed in iRidge et alj (|2006b|). Westend 
encompasses L1448, L1455 and other dark clouds in the South- West part of Perseus. We would like 
to remark that the criteria adopted to identify the sub-regions in the Perseus Molecular Cloud have 
been chosen because they allow us to find (i) the minimum number of regions needed to improve 
the various correlations and (ii) the maximum number of regions with a statistically significant 
number of data points. 

In Figure[3]we present the average ^^CO and ^'^CO spectra for the whole cloud and each region, 
while in Table [2] we show the main properties derived from the average spectra. The central velocity 
and velocity dispersion are computed from the average ^^CO with a Gaussian fit. We can see how 
different the region averages are from each other and the whole cloud. The gradient in central 
velocity across the cloud and the multiple components of the emission are clearly seen. 



5. Analysis 
5.1. Column Density Determination 



To derive the H2 column density, A^(H2), we assume that: (i) all the hydrogen traced by 
the derived extinction is in molecu lar form; (ii) the ra tio between A^(H) and E{B — V) is 5.8 x 

3.1, to obtain 



10^^ cm ^ mag ^ as determined bv lBohlin et al.l ( 19781 ): and (iii) Ry 



Av 



9.4 X 10^° cm ^ mag 



(2) 



The ratio measured by iBohlin et al.l ()1978l ) was calculated only for low extinctions, and it is not 
clear whether extrapolating to higher extinction regions is wise. In addition, it is well kn own that 
Rv can increase up to values close to 4-6 in dense molecular clouds (see e.g. lDraindl2003l ). but we 
assume Ry = 3.1 because it is the average value derived for the Milky Way and therefore our best 
estimation when this quantity has not been measured in the Perseus Cloud. This value of Ry has 
also been used in all previous works, and therefore it facilitates the comparison. 

To estimate ^^CO column densities we assume that the emission is optically thin and in Local 
Thermodynamic Equilibrium (LTE). To estimate the LTE column density, we have to assume an 
excitation temperature, T^x, and optical depth, r. 

In general, the intensity of an emission line, //me; is 



'-line 



{S-h){l-e 



(3) 



where S is the source function and /q is the initial impinging radiation field intensity. The radiation 
temperature is defined as 
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where is the specific intensity and the filling factor is assumed to be unity. Assuming that 



the source function and initial intensity are black-bodies (with 
respectively, then we can write 



Bi,) at Tex and Tf,g = 2.7 K, 



Tr 



Tn 



1 



1 



^To/Tex — \ gTb/Tfcj, _ 
where Tq = hv/k. 

Assuming that the ^^CO (1-0) transition is optically thick, r ^ oo, and that T„ 



(12 



(5) 



CO) is the 



main beam brightness temperature at the peak of ^■^CO, we can derive the excitation temperature 
using eq. ([5]) 



5.5 K 



In (1 + 5.5 K/(r^,,(i2C0) + 0.82 K)) ' 



(6) 



where 5.5 K= hv['^CO)/kB, with i/(i^CO)=115.3 GHz, the frequency of the ^^CO (1-0) line. 



Assuming that the excitation temperature of the ^'^CO (1-0) line is the same as for the ^■^CO (1 



0) hue, the optical depth of ^^CO (1-0) can be derived from eq. 



r(^'^CO) 



In 1 



T™ax(^=^CO)/5.3 K 



l/(e5.3 K/T.. _i)_o.i6, 
where Tm,aa;(^^CO) is the main beam brightness temperature at the peak of ^^CO. 
The formal error for the excitation temperature, a{Tex)-, is 

[Te./(^„^ax(^2CO) + 0.82)]' 



(7) 



(8) 



1 + 5.5/ (T^a,(i2CO)+0.82) 

where ai2 is the error in the ^^CO peak temperature determination. We estimate that ai2 is 
0.15 2~maa; (^^CO) to account for the calibration uncertainty. 



Using the definition of column density (jRohlfs &: Wilsonlll996l ) and expressions ^ and ([7]), we 
derive the ^^CO column density as 



iV(^'^CO) 



1 



-t(13CO) 



3.0 X 10 



1 



-5.3/Te; 



cm 



(9) 



where the Wi}^CO) is the integrated intensity along the line of sight in units of Kkms~^. This 
approximation is accurate to within 1 5% for rf^^C O) < 2 ( Spitzeiill968l ). and always overestimates 
the column density for r(^^CO) > 1 (jSpitzeiill968l ). 

In the determination of the ^^CO column density, we use the derived excitation temperature 
instead of the dust ten iperature becaus e the dust and gas are only coupled at volume densities 
above ~ 10^ cm~^ (e.g. iGoldsmithI l200ll ) . which are typically not traced by ^^CO and ^^CO (1-0) 
data. Moreover, if the volume density of the gas falls below a few times 10^ cm~^ (the critical 
density of the 1-0 transition), the ^^CO lines are expected to be subthermally excited. 
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5.2. Curve of Growth 



Assuming LT E, the photon escape probabihty for a slab as a function of optical depth, /?(t) 
can be written as (|Tielensll2005l ) 



(1 -exp(-2.34r))/4.68r r<7 
l/4r ln(^] r>7 



(10) 



and the Doppler broadening parameter, b, is related to the atomic weight. A, and the gas temper- 
ature, T, by 



^ = 0.1290a/5 kms-i 
V A 



m 



(11) 



where m is the mass of the observed molecule. 



The curve of growth relates the optical depth and the Doppler broadening parameter, b, with 
the integrated intensity, W, by 



where 



W = j Tmb dv = Tjibf{T) 

/(r)=2 r P{f)df 
Jo 



(12) 



(13) 



in the Rayleigh- Jeans regime. In the above expressions Tmb is the main beam brightness temper- 
ature (which is equal to Tr, for a filling factor of unity). 

Assuming that ^^CO emission is optically thin and that the ratio between ^^CO and ^^CO is 
constant in the region, we can write the optical depth as 



where a is the conversion between ^^(^^CO) and ^^CO optical depth. 



(14) 



6. Results 



6.1. Curve of Growth Analysis 



As shown by iLanger et al.l (jl989l ) in B5, ^^CO and ^'^CO integra ted intensities are correlated 
and seem to be well described by the curve of growth (ISpitzeiill978l ). In Figure [5] we show ^^CO 
and ^^CO integrated intensities for Perseus and the individual regions defined in § [H We perform 
a fit of the ""^^CO integrated intensity with a growth curve 



Wi^'^CO) =TRb 



2[3{f)df 



(15) 
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and present the fit results in Table [3] and Figure O 

The fits for the growth curve are very good, considering the simplicity of the model. However, 
it is clear that the correlation is better in the individual regions than in the complex as a whole, 
with the exception of Weste nd, in which t h e fit i s less good. Our results for B5 nicely agree both 
in shape and amplitude with lLanger et al.l (|l989l ) (red solid line in Figure [5]). 

From the fit results we see that for gas at 12 K (intermediate value between average excitation 
temperature and dust temperature) the derived Doppler parameters for B5, IC348, Bl and Westend 
are in reasonable agreement with ay values listed in Table [2j In NGC1333 and the Shell, a larger 
Doppler parameter is expected because the emission comes from multiple components, as seen in 
Figure [H 



6.2. The X factor: using ^^CO to derive A 



V 



The integrated intensity along the line of sight of the CO (1-0) transition, W{ CO), is often 



used to trace the molecular material. The conversion factor X, 

N(R2) =XW{^^CO), 



is derived, and to compare with previous results (e.g. iDame et al.ll200ll ) it is calculated as 

' 9.4 X 10^0 Ay ' 



X 



W^(i2C0) 



(16) 



(17) 



However, the ^^CO emission saturates at Ay ~ 4 mag in every region, as shown in Figure [6j 
Therefore, we perform the estimation of the X factor in two ways: (i) for all the points, and (ii) 
for only those points with Ay < 4 mag, where the column density can still be traced. 

Figure [6] also shows that there is a threshold value of extinction below which no ^^CO emission 
is detected. To take this into account we fit the linear function 



Ay = Ay\2 + ^2 



9.4 X 1020^ 



(18) 



where Ay\2 is the minimum extinction below which there is no ^^CO emission and X2 is the slope 
of the conversion (comparable to the X factor). The linear fit is performed using the bivariate 
correlated errors and intrinsic scatter estimator (BCES), which takes into account e rrors in both 
axes and provides the least biased estimation of the slopes (jAkritas &: Bershadvlll996l ). The results 
of the X factor and the linear fits are presented in Table [H In Figure [6l we show only the results 
for points with Ay < 4 mag: a dotted line for the standard X factor and a dashed line for the 
linear fit. 

The values derived for X using all the points (as previously done in the galactic determinations 
of X) are in agreement with the mean value of (1.8 it 0.3) x 10^*^ cm^^K^^km^^s derived by 
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Dame et al.l (|200ll ) for the Milky Way. However, this fit is not good in the unsaturated regime 
(Ay < 4 mag). Performing a hnear fit to all of the data, including the saturated emission, can 
give unreasonable solutions, such as a negative minimum extinction needed to produce ^■^CO (1-0) 
integrated intensity. 

The linear fit performed to positions with Ay < 4 mag gives the best estimate for the extinction 
in the unsaturated regions but only provides a lower limit extinction estimate for the saturated 
regimes, while the standard X factor provides a poor description of the data in both saturated 
and unsaturated regimes. The histograms of the errors associated with the conversions derived for 
positions with Ay < 4 mag (bottom panels of Figure (H) show that the linear fit (open histogram) 
provides a more unbiased estimate of the extinction than the X factor (filled histogram). This 
improvement in the precision of the extinction estimate goes along a reduction in the errors. The 
width of a Gaussian fitted to the histograms of the regions implies a typical error of 40% and 25% 
for the X factor and the linear fit, respectively. Performing the same analysis on the histograms of 
the whole cloud gives an error of 59% and 38% for the X factor and the linear fit, respectively. 

However, as we can see from Figure [6] the linear relation is not very accurate for Ay > 4 mag. 
Therefore, following the simplest solution of radiative transfer 



we fit it to our data with 



/ = /n 1 



Ty(i2co) = Wo(l- £-'^(^^-^'=12) ) , (19) 



where Wq is the integrated intensity at saturation, ylfci2 is the minimum extinction needed to get 
^2 CO emission, and k is the conversion factor between the amount of extinction and the optical 
depth. We perform an unweighted fit of the non-linear function to the data, which yield solutions 
that better follow the overall shape of the VF(^^CO) as a function of Ay. Unfortunately, the 
best fits for Westend and Perseus are quite poor and should be regarded with caution. The best 
parameters are listed in Table [5] and are shown in Figure [6] as solid curves. We note that the 
threshold extinction shows a large scatter, even in the more accurate non-linear fit. This suggests 
that different environmental conditions are causing the observed scatter (see §[7]). 



6.3. Using ^•^CO to derive H2 column densities 

We expect the ^^CO (1-0) transition to be optically thin at low extinction. If this is the case, 
a simple linear relation between the integrated intensity of ^^CO and Ay should fit the data, 

Ay = WC^CO) Bvi + Awn ■ (20) 

The fit is done using the BCES algorithm for points with Ay < 5 mag because there is some 
saturation in the emission. The results are presented in Table [6] and in Figure [71 Once again, we 
find that the minimum extinction needed to detect ^^CO and the slope of the linear relation changes 
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between the regions, though the variation is smaller than with the ^^CO data. From Figure [7] we 
can clearly see t he effect of saturat ion around 5 mag of extinction in IC348 and Bl, similar to what 
was reported bv lLada et al.l (|1994| ) in the more distant IC 5 146, with cq n ipara ble linear resolution. 
When comparing our results with the linear fit derived by iLada et al.l (|1994| ) we see that the fit 
parameters are quite different, suggesting that the IC 5146 cloud is quite different from Perseus. 
The linear fit for IC 5146 would indicate that for a region without extinction there is molecular 
emission, suggesting that the linear fit has been affected by points with satur ated emission or tha t 
the fit errors are largely underestimated. Comparing our result for B5 with iLanger et al.l (|l989l ). 
we find that the slopes are slightly different, but this can be reconciled by performing th e linear 
fit taking all the points (B13 = 0.37 ± 0.1, Aw^:^ = 1.45 ± 0.04) as done bv iLanger et al.l (|l989l '). 
Ho wever, the fit performed bv iLanger et al.l (|l989l ) was done using the extinction map derived 
by iBachiller &: Cernichard (jl986l ) , that is limited by Ay values below 5 mag of visual extinction 
(due to a lack of detectable background stars). Also, given that their molecular data has a better 
resolution (1.5') than the extinction map (2.5'), they interpolated the ^'^CO to the extinction 
positions instead of smoothing the data to the same resolution. Finally, the threshold extinction 
value differs from previous measurements mainly because it is hard to accurately d efine the zero 



point for extinction maps derived from optical star counting (as already mentioned, ISchnee et al 
2OO5I reported a difference of ~ 0.8 mag). 



Following the column density determination shown in § 15.11 which takes into account the effect 



of optical depth and excitation temperature, we investigate the relation between A^( '^CO) and A 
fitting a linear function to the data. 



V, 



Av = NC'C0)c + Avi3, 



(21) 



where c and Ayis are the parameters of the fit. From the fit we can derive the ratio of abundances 
between H2 and ^^CO as 

H2 



13C0 



9.4 



10^0 c 



(22) 



The result for the fit over the whole cloud and by regions is shown in Figure El and in Table [7] we 
present the best fit parameters. When performing the fit we estimate the error associated with the 
^^CO column density as 15%, due to the uncertainty in the calibration. 

As shown in Figure El this relation presents a larger s catter, which i s also present in individual 



regions (and it is consistent with previous work; see e.g. ICombeslll99ll ). As already pointed out, 
selective photodissociation and/or chemical fractionation can alter the simple linear relation at 
low Ay values, whereas optical depth may start to be large at high Ay (see e.g. the tendency of 
Ni^'-^CO) to flatten out at > 5 mag in B5, IC348 and Bl). 

The ^^CO abundances derived from the fit present significant variations from region to region 
(see Table [7]). We do not find a correlation between the abundance and the threshold Ayi^. This 
suggests that ^'^CO abundance variations are mainly due to different chemical/physical properties 
in the inner regions of the cloud at Ay > Ay 13. 
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Comparing our results with those in Table [71 we see that the abundances agree very well, 
within the errors, w ith previous values r eported in Perseus, taking also into account the ~ 0.8 mag 
difference found by ISchnee et alj (|2005l ) . This difference is produced by the difficulty in defining 
the zero level of extinction, which is harder in the optical star counting method than using NICER 



The e xtinction threshold derived for Perseus is close to the one derived for p-Oph. bv lFrerking et al 



(|l982l ). but they do not have any data at Ay below 2.5 mag and their determination has been done 
only for points with Ay > 4 mag. 

It is important to note that the numbers cited from previous works do not include the 10-20% 
calibration uncertainty that we do include in our results, and therefore our results are more accurate 
than previous ones. 



6.4. Using C^^O to derive H2 column densities 

We fit a linear relation between the integrated intensity of C^^O and Ay, 

Av = VF(Ci80) + Awifi . 



(23) 



The fit results are presented in Tableland in Figure [H We find that, as in the case of ^^CO 
and ^^CO, both parameters (extinction threshold and slope) vary between regions. However, due 
to the fewer points available per region, the errors are larger than for previous fits. 

Despite the small number statistics, we still can see that C^^ O is fairly l i near u p to at least 
10 mag. When comparing our results with the linear fit derived bv iLada et al.l (| 19941 ). we see that 
the fit parameters are quite different. However, they performed the fit over a wider Ay range (up 
to 15 mag), and it is possible that despite their efforts the fit could have been affected by emission 
with higher optical depth. 



L anger et al.l (|l989l ) derived a fit in B5. Unfortunately, we don't have coverage for B5, and 



therefore no direct comparison can be done. Nevertheless, we find that the threshold extinction 
value derived for B5 is systematically lower than the values derived here for other regions (similar to 
what is found in ^^CO) while the slope is in agreement (within the error bars) with the parameters 
derived here. 



6.5. ^^CO (1 — 0) Excitation Temperature vs. Extinction 

Being collisionally excited, the excitation temperature of ^^CO (1-0) is expected to increase 
as we move from the outskirts of the cloud, where the extinction and volume densities may be 
lower than the critical density for the ^^CO (1-0) transition, to the most extinguished and densest 
regions, where ^^CO (1-0) is in LTE and faithfully traces the gas kinetic temperature. 
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The excitation temperature, derived using ^^CO and eq. [H is shown as a function of the visual 
extinction in Figure [TOl In addition we plot the median dust temperature computed in bins of 
extinction as horizontal lines. When all the points are plotted (left panel of Figure [TOj) there is a 
poor correlation between excitation temperature and extinction. This is a direct result of mixing 
very different environments within the cloud in one plot. In fact, from the right panel of Figure \T0\ 
where the excitation temperature and visual extinction are plotted for individual regions, it is clear 
that the scatter is significantly lower and the excitation temperature rises from ~ 5 K at low Ay 
(~ 2 mag) up to a temperature close to the derived dust temperature in positions with Ay > 4. The 
more quiescent regions (Bl and B5) present a smaller dispersion, whereas more active regions (IC348 
and NGC1333) present a larger spread in the excitation temperature, probably because of the larger 
variation of physical conditions along the line of sight and/or the multiple velocity components. The 
region labeled as "Westend" has a very low excitation temperature when compared with the rest of 
the cloud. As shown in the right panel, in individual regions the dust temperature does not change 
more than approximately one degree except in the "Shell" which shows a steady increase with Ay, 
and in NGC1333 where a peak of the dust temperature is present at Ay ~ 8 mag, probably due 
to the internal heating produced by the nearby embedded cluster. On the other hand, the ^^CO 
excitation temperature ranges between 5 and 20 K. 

It is important to note that almost all the points lie below the median dust temperature of the 
region, indicating that ^^CO is t racing gas at volu me densities well below 10^ cm~^, the lower limit 



to have dust and gas coupling (jGoldsmithI l200ll ) . The average excitation temperature for points 



above 4 mag is 13.8 K, while the standard deviation is 2.3 K. We count the number of positions 
where the ^^CO emission is sub-thermal {T^x < 11.5 K) obtaining that it is ~ 60 %. Westend is a 
region with ^^CO excitation temperature always below the dust temperature. This could be due to 
a lower fraction of high density material compared to the other regions. It is interesting that in the 
regions in the North-East part of the cloud (B5, IC348 and Shell) the dust temperature is closer 
to 17 K while in the South- West part (Bl, NGC1333 and Westend) it is closer to 16 K, suggesting 
variations in the ISRF across the Perseus Complex. 



7. Modeling using PDR code 



To relate the observed variations in the ^^CO and ^^CO lines with ch anges in the physical 



properties of the regions, we use the Meudon PDR code (|Le Petit et al.ll2006l '^. This code includes 



most physical effects by explicit calcul ation; in particular it calculates the CO shielding, unlike the 
majority of other codes available (e.g. iRollig et al.l 120071 ). where fitting formulae are used instead. 



We use the abundances derived by iLee et al.l (|l998l ) (see table [9]) for clouds with high metal 
abundances (more appropriate for the material traced by ^^CO and ^^CO) a ^-^C/^^C abundance 



Available through |http : / /axistote . obspm . f r/MIs7] 
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ratio of 80 and a cosmic-ray ionization rate of C = 10~^^ (a value 6 times larger than the 
adopted one has also been considered, but found not to change the results by more than 10%; see 
Appendix). To rep roduce the observ ed [H2/^^C0] ratio we increased the ^^C abundance by a factor 
of 1.8 compared to lLee et al.l (|l998l ). To create curves of integrated intensity as a function of Ay, 
we run PDR models with different extinctions {Av= 0.5, 1, 1.5, 2, 2.5, 3, 4, 5, 8, 10 mag). Each PDR 
calculation is performed over a grid of parameters, assuming a slab of constant density illuminated 
on both sides: (1) The turbulent velocity is fixed to the Doppler parameter derived from the curve 
of growth fit assuming Tr = 12 K for each region (see Tableland § I6.ip . (2) The volume density 
n = 10^, 5 X 10^, 10^, 5 X 10^ cm~^. (3) The radiation field is x = 1 and 3 times the standard Draine's 



radiation field (jDraindll978l ) (other x values have been explored and reported in the Appendix). 



The observed values of Wi^^CO), W{^^CO) and W{^'^CO)/WC^CO) as a function of Ay are 
compared with the models results in Figure [TTl To perform this comparison, the PDR code output, 
I = j lydv, is converted using the definition of brightness temperature (eq. [4]): 



W 



and the ratio can be expressed as 



Tndv 



Ic^ 



(24) 



Ty(i3C0) 



z^(i2C0) ) /(13C0) 



(25) 



First, we note that PDR models are reasonably good in reproducing CO and CO ob- 
servations, if one allows for variations in densities along different lines of sight. One exception is 
NGC1333, where the points with "excess" in ^^CO emission are associated with positions just South 
of the NGC1333 stellar cluster, where a second velocity component is observed. This produces less 
saturation in the ^^CO emission. Secondly, in the case of ^^CO, the Doppler parameter (shown in 
the bottom panels of Figure [TT]) does not affect the results of PDR models within 10%, whereas 
for ^^CO, differences of <40% are present above 2 mag of visual extinction. In fact, the ^^CO 
line is optically thick, thus an increase in the line-width produces a significant increase in the line 
brightness because of the optical depth reduction. The insensitivity of ^^CO integrated intensity 
to variations of the turbulent line-width leaves the density and radiation field as the two possible 
causes of the observed differences between the regions. 

In Figured!] we show only the effects of density with x = 1- Small variations (within a factor 
of 3) of the radiation field intensity slightly shift the curves to higher or lower Ay by 1 mag, if the 
radiation field is larger or smaller, respectively. The main conclusions for individual regions are 
listed below: 



• B5: this is the most quiescent region in the whole Perseus complex. The best fit for the 
VF(^^CO) is a PDR model of a cloud with a narrow range of density change, between 5 x 10'^ 
and 1 X 10^ cm~^, values consistent with previous analysis. A similar density range is found to 
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reproduce W{^^CO) at Ay < 4 mag, whereas 5x 10^ cm~^ is more appropriate for the emission 
at larger extinctions (higher density regions are mainly responsible for the ^^CO emission at 
Ay > 4 mag). We note that the integrated intensities of the two CO isotopologues approach 
(with increasing Ay) the lowest values in the whole sample, suggesting that denser material 
is hidden to view because of photon trapping in the narrow range of velocities observed in 
this region. 

IC348: in this cluster-forming region, the density spread is larger than in B5, with values 
ranging between few x 10^ and ~ 1 x 10^ cm~^. Similarly, the ^^CO at low Av{< 4 mag) 
are reproduced by models with 10'^ — 10^ cm~^, whereas densities larger than 5 x 10^ cm~^ 
are needed above 4 mag. However, the points at Ay > 7 mag and T^(^^CO) > 12 Kkms"^ 
lie next to the embedded cluster, so that local heating and enhanced turbulence are probably 
increasing the ^^CO brightness (the flattening of the W{^^CO) vs. Ay curve is in fact less 
pronounced than in the case of B5, suggesting less optical depth). Thus, proto-stellar activity 
is locally affecting the ^'^CO emission, but not the ^^CO. 

Shell: this region shows the largest spread in density for ^^CO at all values of Ay from 
1 X 10^ to about 3 x 10^ cm ^. At Ay > 5 mag the data groups around two separate values 
of W{^^CO): 35 and 55 K kms"^ Th e forme r group is associated with the (outer) shell 



reported in dust emission by lRidge et al.l (j2006bl ). whereas the latter is located in the inner 



part of the shell, maybe exposed to a larger radiation field causing more dissociation of CO 
molecules. The ^^CO emission appears similar to IC348, with the exception of points located 
at low W(^^CO) (below 4 Kkms~^) and Ay > 4 mag which are again associated with the 
inner part of the Shell. This also suggests some further destruction process for the ^^CO. 

• Bl: the ^^CO emission is consistent with material at densities between 5 x 10'^ — 3 x 10^ cm^^. 
Compared with B5, the Bl region shows brighter ^^CO lines at lower (as well as higher) 
extinction. This is probably related to the larger Doppler parameter of Bl. Several data 
points at low Ay (<2 mag) are well reproduced by PDR models of dense clouds. The need 
of high densities at low Ay also appears in the ^^CO panel. 

• NGC1333: this is the most active star forming site in the whole Perseus cloud and the be- 
haviour of the ^^CO and ^'^CO integrated intensities as a function of Ay is in fact significantly 
different when compared to the other regions. First of all, the densities required in the PDR 
code to match the data are mostly above 10^ cm~^, both for the ^^CO and ^^CO emission. 
Secondly, the saturation of the ^^CO line becomes evident only at Ay > 6 mag (unlike 
~ 4 mag, as in the other regions). Here, similarly to what is seen in Bl, non-thermal motions 
driven by the embedded protostellar cluster are broadening the CO lines, allowing photons 
from deeper in the cloud to escape. The effect is more pronounced than in Bl, consistent with 
the fact that NGC 1333 has the largest Doppler parameter among the six regions. We further 
note that, unlike in IC348, the ^^CO (1-0) integrated intensity is also affected by the internal 
star formation activity, significantly reducing the saturation and enhancing the brightness at 
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large Ay- Internal motions, likely driven by protostellar outflows, are thus more pronounced 
in NGC 1333 than in IC348, likely because of the larger star formation activity. 

• Westend: this is the only region where no data points are present at Ay > 6 mag, and the 
^^CO, as well as the ^^CO, integrated intensity shows a large scatter between Ay of 1 and 
6 mag. These two facts are consistent with an overall lower density and probably clumpy 
medium, where relatively small high density clumps are located along some lines of sight, 
whereas a significant fraction of the data (19% of points in ^-^CO) can be reproduced by 
uniform PDR model clouds with densities below 5 x 10^ cm~'^. 



In general, model predictions for ^^CO (1-0) can only reproduce well the observed emission at 
low extinction {Ay < 3 mag). The complex structure of active star forming regions, in particular 
density and temperature gradients as well as dumpiness along the line of sight (all phenomena not 
included in the PDR code) can of course contribute to the deviations from the uniform PDR models. 
We point out again that the largest Doppler parameters, i.e. the largest amount of non-thermal 
(turbulent?) motions, are present in active regions of star formation, so their nature appears to 
be linked to the current star formation activity and not to be part of the initial conditions in the 
process of star formation. 

In the bottom panel of Fig. dH the W{^^CO)/W(^^CO) ratio is shown as a function of Ay 
for the six regions. The PDR models appear to reproduce well the integrated intensity ratio, for a 
broad range of Ay. As we just saw, the ^^CO data preferentially trace higher density material than 
^■^CO lines, so the black squares show the ratio between the ^^CO and ^^CO emission as predicted 
by PDR models with 5 x lO^cm"^ for ^^CO and 1 x lO^cm"^ for ^"^CO lines. One thing to note in 
these plots is the large fraction of points at low Ay and low W(^'^CO)/W{^^CO) which lie below 
the PDR model curves, in particular for B5, NGC 1333, and Westend, but they lie above the 
black squares showing that the PDR model with different densities can reproduce all the emission. 
Another way to reproduce these data points is by decreasing the interstellar radiation field by a 
factor of a few. Alternatively, it is possible that these lines of sight intercept material where the 
^^C carbon is stih pa rtially in ioniz ed form, so that the reaction ^^C+ + ^^CO ^^C^ + ^^CO + AE 



(with AE/k = 37 K: IWatsonI 119771 ) can proceed and enhance the CO abundance relative to CO. 



We finally note here that iBell et al.l (|2006l l have theoretically investigated the variation of the 
X factor using UCL_PDR and Meudon PDR codes. They argue that variations in X can be due to 
variations in physical parameters, such as the gas density, the radiation field and the turbulence, 
in agreement with our findings. 



8. Summary and Conclusions 

Using the FCRAO ^^CO, ^^CO and C^^O data, and a NICER extinction map produced by 
COMPLETE we perform a calibration of the column density estimation using ^^CO, ^^CO and 
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emission in Perseus. We report the following results: 

• We find a parameter space, VLSRC^CO)-Td, in which different spatial regions of the Perseus 
Molecular Cloud Complex also cluster, and we designate six regions within the complex (see 
Figures [2] and [3]) . We note that the dust temperature decreases by about 1 K from North 
East to South West. 

• The ^^CO data can be modeled with a curve of growth. The fit parameters vary between 
the six regions and this causes much of the scatter in the I^(^^CO) vs. W{^^CO) plot of the 
whole Perseus Complex. The parameters derived from the fits agree with a previous study of 
a sub-region of Perseus (B5) to within errors (see Figure [5]). 

• The X factor, X = N(R2)/W{^'^C0), is derived from linear fits to the data both for the 
whole Perseus Complex and for the six regions. The ^■^CO saturates at different intensities 
in each region, depending on the velocity structure of the emission, the volume density and 
radiation field. When the linear fit is done only for the unsaturated emission, the X factor 
is smaller than that derived for the Milky Way. However, larger values are obtained (closer 
to that found in the Milky Way) if all the points are included in the fit (see Figure [6] and 
Table H]). The most active star forming region in Perseus (NGC 1333) has the lowest X factor 
and the largest ^^CO abundance among the six regions. 

• The gas excitation temperature varies from 4 K to 20 K, it increases with Ay, and it is 
typically below the dust temperature at all Ay- This can be explained if a fraction (~ 60%) 
of the ^^CO (1-0) lines is sub-thermally excited, i.e. if the ^^CO- emitting gas has volume 
densities below ~ 3 x 10^ cm~^. 

• The column density of ^^CO is derived taking into account the effect of optical depth and 
excitation temperature. We find that the threshold extinction above which ^^CO (1-0) is 
detected is larger than has previously been reported. However, the fractional abundances 
(w.r.t. H2 molecules) are in agreement with previous determinations. The difference with 
previous works is due to the superior zero-point calibration and larger dynamic range of the 
NICER extinction map, as compared to those derived from optical star counting (see Figure [8] 
and Table [7]) . 

• ^^CO abundance variations between the regions do not correlate with the extinction threshold 
Avi3, suggesting that the main cause of the variation is likely due to the chemical/physical 
properties of shielded molecular material deeper into the cloud. The ^■^CO (1-0) and ^'^CO (1- 
0) lines saturates at Ay > 4, 5 mag, respectively, whereas C"*^^© (1-0) line do not show signs 
of saturation up to the largest Ay probed by our data (10 mag). 

• Using the Meudon PDR code we find that the observed variations among the different regions 
can be explained with variations in physical parameters, in particular the volume density 
and internal motions. Large Doppler parameters imply large values of the CO integrated 
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intensities (as expected for very optically thick lines) and are typically found in active star 
forming regions (the largest values of the Doppler parameter and Vl^(^^CO) being associated 
with NGC 1333, the most active site of star formation in Perseus). On the other hand, 
quiescent regions such as B5 appears less bright in CO and only show a narrow range of CO 
integrated intensities as a function of Ay. This is likely due to the fact that the photons 
emitted from the higher density regions located deep into the cloud have similar velocities 
relative to the outer cloud envelope traced by ^^CO, so that they are more easily absorbed. 
Thus, turbulent (or, more generally, non-thermal) motions appear to be a by-product of star 
formation, more than part of the initial conditions in the star formation process. 

This work has shown that local variations in physical conditions significantly affect the relation 
between CO-isotopologue emission and Ay, contributing to the observed scatter. The use of a 
standard X factor, 1.8 x lO^'' cm~^ K^^ km^^ s, produces an overestimation of the cloud's mass by 
~ 45% when compared to the mass derived from the extinction map, while the lower limit for the 
mass derived using the linear fit to the unsaturated points underestimates the mass by a ~ 15%. 
The X factor (as well as the ^^CO fractional abundance) depends on the star formation activity, 
with lower values associated with the more active (and turbulent) regions. Extinctions measured by 
using ^^CO and previous conversions from the literature are typically underestimated by ~ 0.8 mag, 
so that more shielding is needed to produce the observed ^^CO compared to previous findings. 

JEP is supported by the National Science Foundation through grant t^AF002 from the Asso- 
ciation of Universities for Research in Astronomy, Inc., under NSF cooperative agreement AST- 
9613615 and by Fundacion Andes under project No. C-13442. This material is based upon work 
supported by the National Science Foundation under Grant No. AST-0407172. PC acknowledges 
support by the Italian Ministry of Research and University within a PRIN project. 



A. Effect of X and ( variation 

In Sect. Owe have explored how changes in volume density and Doppler parameter affect the 
i^CO (1-0) and ^^CO (1-0) integrated intensities (Wi^^CO), W{^^CO)) predicted by the Meudon 
PDR code. Here we show the effects of variations in the interstellar radiation field intensity, in units 
of Draine's field (x), and the cosmic-ray ionization rate {() on VF(^^CO) and VF(^^CO). Fig. [T2] 
shows the results of this parameter space exploration in the particular case of B5 and volume 
density of 5 x 10'^ cm~'^ (similar results apply to the other regions and different densities). The 
upper panels display the model results for x values of 0.5, 1, 5 and 10: the main change is visible 
at Ay < 4 mag, with a shift of the threshold extinction for ^^CO and ^'^CO emission from about 1 
to 3 mag for an increase of x from 0.5 to 10, respectively. 

The cosmic-ray ionization rate used in the PDR models described in Sect.[7|is = 1 x 10^^^ s^^. 
This value is quite uncertain and I Dal gar nd ( 20061 ) suggests a higher rate of 6 x 10^^^ s^^ for 
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molecular clouds (see also Ivan der Tak fc van Dishoeckll200d ). In the bottom panels of Fig. \T2\ we 
show the predicted W{^'^CO) and VF(^^CO) curves for C = 6 x 10"^^ s~^. The larger C, value does 
not affect the CO emission at Ay < 3 mag, and only changes the integrated intensities by about 
30%. Thus, the effect is not large enough to explain the largest W{^^CO) values. 
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Table 1. Clusters Regions Removed 



Region 


Center R.A. 


Center Decl. 


Box Size 




(deg) 


(deg) 


(deg2) 


IC 348 


56.088 


32.171 


0.480x0.426 


NGC 1333 


52.212 


31.483 


0.397x0.496 



Table 2. Typical Properties of the Regions from Average Spectra 



Region 


H/(i2cO) 


Trnax i^"^ CO) 






iy("co) 


Tmax {^'^ CO) 


Visfl("CO) 






t(13CO) 




(K km s-i) 


(K) 


(km s"-*-) 


(km s^-"-) 


(K km s-i) 


(K) 


(km s^^) 


(km s^-') 


(K) 




B5 


8 


3.2 


9.84 


0.90 


1.6 


0.99 


9.99 


0.64 


12 


0.31 


IC348 


9 


3.0 


9.01 


1.18 


2.2 


0.93 


8.99 


0.95 


12 


0.35 


Shell 


10 


2.6 


8.73 


1.49 


2.0 


0.58 


8.72 


1.30 


11 


0.28 


Bl 


12 


2.6 


6.71 


1.55 


1.6 


0.99 


6.83 


0.99 


11 


0.33 


NGC1333 


15 


3.0 


6.68 


1.85 


2.6 


0.73 


7.06 


1.38 


11 


0.31 


Westend 


10 


2.2 


4.20 


1.96 


2.0 


0.65 


4.58 


1.19 


9 


0.44 


Perseus 


11 


1.7 


7.26 


2.65 


2.2 


0.39 


7.67 


2.15 


11 


0.34 



Note. — W = integrated intensity. Tmax = peak brightness temperature. Visfl = centroid velocity, cry = velocity dispersion (= Au/ V8 In 2, where 
Av is the full width at half maximum). T^x = excitation temperature derived from ^-^CO. t = ^^CO optical depth. The l-cr uncertainty for integrated 
intensity, peak brightness and excitation temperature is estimated between 15 and 30%. 
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Table 3. Parameters for Growth Curve Fit 



Region 


a 








(Kkms-i) 


B5 


0.61±0.03 


20.5±0.5 


IC348 


0.246±0.007 


33.5±0.6 


Shell 


0.223±0.009 


40±1 


Bl 


0.33±0.01 


36.8±0.7 


NGC1333 


0.130±0.006 


67±2 


Westend 


0.44±0.03 


24.7±0.9 


Perseus 


0.260±0.004 


35.9±0.4 
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Table 4. Linear fits to ^^CO 



Region X/IO^" Avi2 Xa/lO^" 

(cm~^ km"-' s) (mag) (cm^^ K"-' km"-' s) 



Fit performed to the whole dataset 



B5 


2±1 


0.0±0.2 


1.5±0.1 


IC348 


3±2 


0.51±0.08 


1.61±0.07 


Shell 


2±1 


l.OiO.l 


1.09±0.08 


Bl 


1.4±0.8 


-0.6±0.1 


1.38±0.09 


NGC1333 


0.9±0.3 


-0.5±0.1 


0.93±0.07 


Westend 


1.2±0.5 


-1.8±0.3 


1.8±0.4 


Perseus 


2±1 


-0.33±0.07 


1.38±0.05 




Fit performed to points where Av 


< 4 


B5 


2±1 


0.97±0.06 


0.88±0.09 


IC348 


3±2 


1.60±0.04 


0.76±0.06 


Shell 


2±2 


1.5±0.1 


0.8±0.1 


Bl 


1.4±0.9 


1.18±0.04 


0.56±0.04 


NGC1333 


0.9±0.3 


0.60±0.06 


0.57±0.06 


Westend 


1.1±0.4 


-0.4±0.2 


1.2±0.3 


Perseus 


2±1 


0.92±0.04 


0.72±0.04 


Previous -works 


Galaxy (1) 


2.8±0.4 






P-Oph (2) 


1.8±0.1 






Galaxy (3) 


1.8±0.3 






Pipe (4) 




2.02±0.02 


1.06±0.02 



Ref erences. — (1) iBloemen et"al] (|l986h : (3) IPrerking et"ZI (|l982l '): (3) IPame et al] (|200ll ): (4) iLombardi eTal 
( 200i ): 
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Table 5. Parameters for CO Fit 



Region 




k 


Akl2 




(Kkms-i) 


(mag^"'^) 


(mag) 


B5 


30.92 


0.553 


1.063 


IC348 


39.27 


0.350 


1.374 


Shell 


73.31 


0.139 


0.851 


Bl 


43.08 


0.691 


1.309 


NGC1333 


67.12 


0.374 


0.748 


Westend 


30.45 


0.529 


-0.160 


Perseus 


42.29 


0.367 


0.580 



Table 6. Results of Linear Fit to W{^^CO) 



Region 




Biz 


Reference 




(mag) 


(magK~^ km^"*^ s) 




B5 


1.56±0.03 


0.323±0.009 


1 


IC348 


1.99±0.03 


0.36±0.01 


1 


Shell 


1.90±0.06 


0.40±0.02 


1 


Bl 


1.64±0.03 


0.296±0.009 


1 


NGC1333 


1.19±0.03 


0.26±0.01 


1 


Westend 


0.75±0.06 


0.44±0.02 


1 


Perseus 


1.46±0.02 


0.345±0.006 


1 


Previous works 


IC 5146 


-2.6±0.3 


1.4±0.1 


2 


B5 


0.54±0.13 


0.39±0.02 


3 



References. — - (1) T his work, linear fit for Ay < 5; (2) iLada et al.l (119941 ). linear fit for Ay < 5; 
(3) banger et all (Il989l ) 
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Table 7. Results of Linear Fit to N{^^CO) 



Region 


^V13 


c 


[H2/13CO] 


Reference 




(mag) 


(mag cm~^) 


xlO^ 




B5 


1.62±0.04 


(4.0±0.2)xl0"^^ 


3.8±0.2 


l,a 


IC348 


2.15±0.04 


(4.4±0.1)xl0~^^ 


4.1±0.1 


l,a 


Shell 


2.23±0.05 


(4.3±0.1)xl0-^'^ 


4.1±0.1 


l,a 


Bl 


1.68±0.04 


(4.0±0.1)xl0-^6 


3.8±0.1 


l,a 


NGC1333 


1.44±0.04 


(3.0±0.1)x 10-1*5 


2.8±0.1 


l,a 


Westend 


1.19±0.05 


(5.2±0.2)x 10-1*5 


4.9±0.2 


l,a 


Perseus 


1.67±0.02 


(4.24±0.07)x 10-1*5 


3.98±0.07 


l,a 


Previous works 


Perseus B5 


0.5±0.1 


(3.8±0.2)x 10-1*5 


3.6±0.2 


2,b 


Perseus 


0.8±0.4 


(4.0±0.8)x 10-1*5 


3.8±0.8 


3,b 


L1495 


0.3±0.3 


(4.5±0.6)xlO"i*5 


4.2±0.6 


4,b 


L1517 


0.3±0.5 


(6±2)x 10-1*5 


5±1 


4,b 


p-Oph 


1.6±0.3 


(3.7±0.4)x 10-1*5 


3.5±0.4 


5,c 


Taurus 


1.0±0.2 


(7.1±0.7)x 10-1*5 


6.7±0.7 


5,c 



Note. — (a) Ay derived from NIR colors, (b) Ay derived from star counting; (c) Ay derived 
using spectra. 



Duvert et al.l (|l98m : (5) iFrerking et al.l mm 



References . — ( 1) T his work; (2) iLanger et al.1 \l989i ): (3) iBachiller k Cernicharol (11983); (4) 
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Table 8. Results of Linear Fit to W{C^^O) 



Region 




-Bis 


Reference 




(mag) 


(mag km^"*^ s) 




IC348 


2.1±0.4 


5±11 


1 


Shell 


2.7±0.3 


4±7 


1 


Bl 


2.1±0.1 


2.5±0.9 


1 


NGC1333 


1.7±0.2 


3±2 


1 


Westend 


1.6±0.3 


3±5 


1 


Perseus 


2.4±0.1 


2.9±0.9 


1 


Previous works 


IC 5146 


-0.7±0.3 


lOibl 


2 


B5 


1.40±0.22 


1.8±0.13 


3 



Ref erences. — (1) This work; (2) ILada et alJ (|l994l ). linear fit for Ay < 15; (3) iLanger et al 



(119891) 



Table 9. Initial chemical abundances with respect to total Hydrogen 



Element Abundance 



He 


0.1 




C+ 


1.307 X 10 


-4 




1.633 X 10 


-6 


N 


2.14 X 10" 


5 





1.76 X 10- 


4 


s+ 


8.00 X 10- 


6 


Fe+ 


3.00 X 10- 


7 




Fig. 1. — Error in the Perseus extinction map derived with NICER at each pixel versus the corre- 
sponding Ay- Only the pixels used in this study are displayed here. 
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Fig. 2. — Central velocity compared with the dust temperature. Perseus is divided into six regions 
that are defined mainly in the space of physical parameters. This separation allows a better 
understanding of the cloud. 
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Right Ascension {J2000) 



20"0.00° 



Fig. 3. — (Top) The extinction map derived using NICER. The green border is the region ob- 
served in ^^C O and ^^CO (1-0) by COMPL ETE, while the red border is the region observed in 
C^^O (1-0) by iHatchell k van der TakI (|2003l ). {Bottom) The re-gridded molecular data that fulfil 
the requirements detailed in §[3] are shown as boxes. Each of the defined regions are presented in a 
different color. B5 is green, IC 348 is magenta, the Shell is cyan, Bl is orange, NGC 1333 is blue 
and Westend is red. 
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Fig. 4. — Average ^^CO (1-0) and ^^CO (1-0) spectra for all Perseus and the six sub-regions. 
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Fig. 5. — Integrated intensity of ^■^CO is plotted against the integrated intensity of ^^CO. The left 
panel shows all data used while right panel shows each region separately, using the same colors as 
in Figure [3l The median of the 1-sigma errors are shown in bott om right of each plo t. Solid lines 
are the growth curve fit, while the red curve in B5 is the fit from lLanger et al.l (|l989l ). 
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Fig. 6. — Top: Integrated intensity of ^^CO and visual extinction derived using NICER. The left 
panel shows all data used and the right panel shows each region separately, using the same colors 
as in Figure [3l Solid lines show the best fit of eq. (fT9j) , dotted and dashed lines are the standard 
W(^'^CO) X = iV(H2) and straight line fits are only for points below Av=^, respectively. 5ottom: 
Histograms of the difference between the observed and expected Ay for all data displayed in top 
panels. Filled and open histograms are for the standard X factor and linear conversion, shown as 
dotted and dashed lines in upper panels, respectively. 
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Fig. 7. — Left panel shows CO Integrated intensity compared with the extinction. The median 
of 1-sigma errors are shown in the bottom right corner of each plot. The best linear fit is shown 
with a solid line. The right panel shows same plots as the left panel, but for the individual regions. 




Fig. 8. — The left panel shows CO column densities as a function of visual extinction. The 
median of the 1-sigma errors are shown in bottom right of each plot. The excitation temperature 
is estimated using the ^^CO peak temperature, assuming this transition to be optically thick. The 
best fit is shown with a solid line. The right panel show the same as the left panel, but separated 
into different regions. 
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Fig. 9. — The left panel shows C O integrated intensities as a function of visual extinction. The 
median of 1-sigma errors are shown in the bottom right corner of each plot. The best linear fit 
is shown with a solid line. The right panel shows the same as the left panel, but separated into 
different regions. 




^ 10 



. ' 1 ' 1 ' 1 ' 1 ' 1 ' . 


. ' 1 ' 1 ' i_ ' 1 ' 1 ' 


. ' 1 ' 1 ' 1 ' 1 ' 1 ' . 


r.! 


mP \ - 

■' ica« 


: ' 1 ' 1 ' 1 ' 1 ' 1 ' : 


: ' 1 ' 1 ' 1 ' 1 ' 1 ' 

■ » 


: ' 1 ' 1 ' 1 ' 1 ' 1 ' : 


Ei : 

■ . 1 . 1 . 1 . 1 , 1 . ■ 


r,: 

■ . 1 . 1 . 1 . 1 . 1 . 


r^" + 

■ . 1 . 1 . 1 . 1 . 1 . ■ 



C 2 4 G a 10 



Fig. 10. — Excitation temperature derived using ^^CO plotted against the visual extinction at the 
same position. The median of the 1-sigma errors are shown in bottom right of each plot. The solid 
line is the median dust temperature for each sample within the extinction bin. The left panel shows 
all data used while the right panel shows them separated by regions. 
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Fig. 11. — Results of the PDR modeling. Upper, middle and bottom panels show the comparison 
between W{^^CO), W{^^CO) and the ratio W{^^CO)/WC^CO) with Ay, respectively. Black and 
grey lines are derived from PDR models. Unconnected black squares show the ratio between ^^CO 
and ^'^CO emission predicted by the PDR model with densities of 5 x 10'^ and 1 x 10^ cm~^, 
respectively. 
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Fig. 12. — PDR results for the particular case of B5 and different values of x C (iii units of 
10^^'' s"^). Black solid, black dotted, gray solid and gray dotted lines show the PDR model results 
for X = 0.5,1.0,5.0 and 10 times the standard Draine's radiation field, respectively. Variations in 
X mainly affect the low Ay emission, whereas different C values slightly change the emission values 
at Ay > 4 mag. 



